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Abstract 

We revisit the bremsstrahlung process of a superluminal neutrino motivated by 
OPERA results. From a careful analysis of the plane wave solutions of the superlu- 
minal neutrino, we find that the squared matrix elements contain additional terms 
from Lorentz violation due to the modified spin sum for the neutrino. We point out 
that the coefficients of the decay rate and the energy loss rate significantly depend on 
the details of the model, although the results are parametrically similar to the ones 
obtained by Cohen and Glashow [1]. We illustrate this from the modified neutral 
current interaction of neutrino with Lorentz violation of the same order as in the 
modified dispersion relation. 
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1 Introduction 



Lorentz invariance is one of the cornerstones of the modern quantum field theory, and 
it was completely compatible with all previous experiments and observations. Recently, 
an intriguing result has been presented by the OPERA collaboration [2j, claiming that 
the muon neutrino speed exceeds the speed of light by (f — c)/c ~ 2.37 ± 0.32^q;24 x 
10~^. This result has been recently confirmed by a test performed using short-bunch wide- 
spacing beam in the revised version of Ref . [2] . There are similar measurements from other 
oscillation experiments, for example, MINOS, resulting in {v — c)/c = 5.1 ± 2.9 x 10~^ |3], 
but the OPERA is the first experiment which observed a positive {v — c) for neutrinos 
with high statistical significance of about 6.2cr. There is no statistically significant energy 
dependence of superluminality at OPERA where the results for the low- and high-energy 
samples with averaged neutrino energies, 13.8 GeV and 40.7 GeV, respectively, are in 
agreement [2]. 

Such a drastic result is seemingly in contradiction with a set of other neutrino obser- 
vations 0. The detection of neutrinos emitted from the SN1987A supernova [SI El E] puts 
a stringent bound on the electron anti-neutrino speed, |f — c|/c < 2 x 10^^. Furthermore, 
the observation of neutrino oscillations demands that neutrino velocities for different neu- 
trino flavors should be equal up to \vi — Vj\ < 10~^^ for i ^ j, otherwise the coherence is 
lost and the oscillation pattern is smeared [H E]. This contradiction could be explained 
in two ways. One is to make the neutrino velocity energy dependent [lOl Hlj, because the 
typical neutrino energy is 10 MeV for the supernova neutrinos, and is about 28 GeV for 
the CNGS neutrino beam, used by OPERA. Another is to make the neutrino superluminal 
only within the Earth radius [121 El 1131 [IH [15] , or only inside matter [16] . A different 
route to the solution could be taken by considering models with energy non-conservation 
or deformed Lorentz invariance [T7] in the neutrino interactions. But this would be hard 
to formulate in the language of ordinary quantum field theory so we do not pursue this 
option in this article. 

In any model explaining the OPERA results it is important to check, whether the 
creation, propagation, and detection in the OPERA setup can be explained. There are 
very strong statements that invalidate most of the proposed models for the OPERA results. 
First, the superluminal neutrino can radiate electron-positron pairs (in a way analogous 
to the Cherenkov radiation) [1] , thus loosing energy before reaching the detector. Second, 
the decay of a fast moving pion is modified, and even the initial neutrino spectrum should 
have a strong cutoff at energy, which is below the average energy detected by OPERA 
[TH| dSl I2D] • Both of these results rely on the following assumption — the only thing modified 
in the theory is the dispersion relation of the neutrino. As the neutrino speed is given by 
the derivative of the dispersion relation v = dE/dp, a constant neutrino speed at OPERA 
means that the neutrino dispersion relation has the form E = vp = {1 + S/2)p with 
(5 ~ 5 X 10~^. A stronger claim, based on the result of Ref. [Ij was made by the analysis 
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of the ICARUS results [21]. The ICARUS detector, while not being able to measure the 
arrival time of the neutrinos, can carefully measure their energy spectrum. The comparison 
of the expected and measured spectra provides strict bounds on the neutrino speed, because 
the energy loss by electron-positron radiation pjj would significantly change the spectrum. 

In this article we re-analyze the bremsstrahlung process of a superluminal neutrino 
that has been done in Ref. [1]. The main point of the previous analysis follows from the 
kinematical possibility for the neutrino with non-standard dispersion relation to "decay" 
into other particles. For example, for the process v — )■ z/ + + e~, the "masses" (or 
squares of the four momenta) of the initial and final neutrinos are different, enabling the 
process to take place. However, the exact calculation for the process rate is more involved. 
Specifically, when the dispersion relation for neutrino is modified, one has to use the 
modified plane wave solutions for the neutrino. In turn, the spin sum for the final neutrino 
gets modified, leading to the additional terms of the order of 5 in the expression for the 
squared matrix elements, as compared to the previous calculations [I]. Moreover, when 
there are modifications of the same order in the electroweak interaction vertex of neutrino 
as in the neutrino dispersion relation, there appear more terms of the similar order in the 
squared matrix elements too. All these effects are added up to give a nontrivial result, 
which depends on the additional modifications of the same order. This is due to the fact 
that the squared matrix elements (obtained by the standard rules for the spin sum for 
Lorentz invariant fermions) is only of the order of 5^ in the kinematically allowed region. 
As a result, due to various cancellations in the matrix element, the final probability of the 
bremsstrahlung process depends on the details of the Lorentz violation in the model. 

We explicitly construct (at the level of Fermi four fermion interactions) two models with 
broken Lorentz symmetry. Both models have a common property that the neutrino kinetic 
term contains a Lorentz violating term as inspired by the modification of the metric for the 
neutrino [12] . The difference is that one of the models keeps the interaction terms Lorentz 
invariant while the other model introduces a similar Lorentz violation in the electroweak 
neutral current of neutrino too. 

In section [2] we introduce the Lagrangians for the models. Then, in section |3] we obtain 
the free solutions for the neutrinos and present the rules for "summation" over the spin 
states. Consequently, in section H] we provide the detailed calculation of the neutrino decay 
width and the rate of energy loss. Finally, conclusions are drawn. 



2 Models 

We first define the framework for the calculation of the decay (bremsstrahlung) of a super- 
luminal neutrino into an electron-positron pair. Following the ideas in [SI [22] the lowest 
order Lorentz violating operator in the Lagrangian for the massless (Weyl) fermion looks 
like 

li = iv^^~g,,d\\-^'')v, (1) 
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where the Lorentz violating "metric" can be chosen as 



g^^ = diag(l, -v, -v, -v) 



(2) 



with the neutrino speed v = 1 + 6/2. Quite obviously, this action gives rise to the super- 
luminal neutrino dispersion relation, E = v\p\. In the limit of Lorentz invariance, we get 
g^,^ = diag(l, -1, -1, -1) = g^^. 

The next component of the model is the interaction term. We will not go in the details 
of how the model emerges from the underlying electroweak gauge theory, because this 
would lead to complications due to the different velocities for the left-handed electron and 
neutrino of the same multiplet. We will take the purely phenomeno logical approach and 
analyze two types of the four fermion neutral current interaction. The first one will be the 
usual Lorentz invariant one (model I) 



intl 



G p 

7! 



z/7^(l - 75)z/ eY{ve - ae7^)e 



(3) 



The second one is inspired by a "gauge invariant" Lagrangian, where the covariant deriva- 
tive enters in the same way as in ([1]) (model II) 



-^int2 — — ;= 
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(4) 



In fact, this interaction term does not follow from the gauge invariant SM, as far as elec- 
troweak gauge symmetry is unbroken. In this work, we assume this possibility but do not 
consider a microscopic model with Lorentz symmetry /electroweak symmetry breaking for 
that. 



3 Free solutions for the neutrinos and spin sums 

The action ([T]) leads to the Dirac equation of the form (in momentum representation and 
two component form for simplicity) 

(i^ffO - vfa')x = 0, (5) 

where cr^, a* are the unit 2x2 matrix and Pauli matrices, respectively. Let us also (using 
invariance under 0(3) spatial rotations) align the momentum along the 3rd spatial axis. 
Then we immediately get two solutions 

E = vp^ with X = 72E , (6) 

E = -vp^ with X = V2E , (7) 

where \/2E is the standard overall normalization. The first solution corresponds to the 
neutrino and the second solution to the antineutrino, both with velocity v. Notice that 
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the dispersion relation is modified, while the spinors have the usual form. The spin "sum" 
for the neutrino (which is in this case trivial, as far as there is only one spin state for the 
neutrino) is 

J2 XsXl = 2E [J = p'^~g,,a^ ^ jfg.^a^. (8) 

s=l/2 ^ ^ 

The important observation here is that the spin sum is not given by the standard expression, 
but with the momentum contracted with the sigma matrices using the superluminal metric. 
The difference is nontrivial in the first order in 5, which is essential for the calculation of 
the decay width. 

The generalization to the usual four component spinors is obvious, 

^ z/,z/, = p^g^^Y = ^, (9) 

s 

where the momentum with the tilde is a shorthand for p^ = g^vP'^ ■ Note, that raising 
and lowering the indices is always done with the normal metric, while the tilde means the 
additional factor of v for the spatial part in the scalar product. 

Now we are ready to evaluate the matrix elements and the squared of them. 



4 Decay width calculation 



Following the standard rules for the decay process, we get for the bremsstrahlung process 
i^/x — ^ ^'/j + e"*" + e~ with four momenta p, p\ qi, and q2, 



i2nl 
2E, 



n 



d^Pf 



(27r)32E; 



(10) 



spin 



Here, note that we do not average over the spin states of the neutrino, which was done in 
the calculation of [1]. In the current setup the neutrino has explicitly only one spin state, 
instead of two spin states for a massive neutrino. The square of the matrix elements is the 
following: in the model I ([3]), 



|2 = ^W^Er. 



and in the model II 



J2\Mn\' = ^M-^g^,g,sE^'. 
Here the individual traces are 

M"'^ = tr 



Eai3 = tr 



^7^(1-75)^(1 + 75)7'' 



(11) 

(12) 

(13) 
(14) 
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There are two differences from the standard calculation with Lorentz invariance: the spatial 
parts of p and p' momenta in M"'^ are multiplied with v, and the indices between M and 
E are contracted with superluminal metric for the model II. 

The results of the multiplication are (in the approximation of a purely axial electron 
neutral current, Oe = —1/2 and fg = 0) 



{pqi){p'q2) + {p'qi){pq2) 



(15) 



and 



{pqi){p'q2) + {p'qi){pq2) 

-{pp'){qiq2) - {pp'){qiq2) 



{pp){qiq2) 



(16) 



Here the tilde always means one factor of v in front of the spatial product, i.e. (pq) = 
pOgO _ . (^pq'j ^ pOqO _ . ^ {pp') = p^p'^ — w^p ■ p', etc. Note that in model II the 
second line in f|T6|) does not vanish. For comparison to our model I, in Ref. [23], the second 
term in the squared amplitude (1151) was missing, the spin average for the initial neutrino 
was taken, and the modified plane-wave solutions for neutrino were not taken into account 
for the spin sum of the neutrino. 

The rest of the calculation is rather straightforward, and consists just of careful inte- 
gration over the final momenta. The safest way is to perform the calculation in the lab 
frame directly. We will only sketch the derivation here. 

First, we perform integration over the momenta of the electron and positron, in the 
limit of zero electron mass (this is fine as far as we are interested in the decay of high 
energy neutrinos) 



qif^q2u-^ [k 

-C/l -C/2 



qi - 92) 



n 2 



(17) 



where k = qi + q2 is the momentum of the electron-positron pair. It is convenient to rewrite 
the remaining integration as the integration over the modulus of the final neutrino momen- 
tum Ip'I and the cosine of the angle 6 between p and p' (in the lab frame). Calculating 
all the scalar products together with the dispersion relations, 
obtain the decay rate as follows. 



f IpI and p = v\p 



we 



gg7]-3^2|p| 



|p'|(i|p'|(icos^^ X 



where I = Xm for models I and II are given by 



Xi = ((t;2 - l)(|p|2 + |pf) - 2\p\\p'\{v^ - cos0))?;'|p||p'|(l - cos^) 

- 2v^(^{v - l)|pp - |p||p'|(t; - cos^)) (^{v - l)|pf - IpIIp'Ki; - cos^)) (19) 
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and 



X, 



II 



t;2|p||p'|(l-cose) 



;i+3t;2)(^(t;2_i)(|p|2+|p'|2)_2|p||p'|(^2_^Qg^^j^2w2|p||p'|(i_cos 
-2(^(t;2-i)(|p|2+|p'|2)_2|p||p'|(^2_^^g^)j^2|p||p/|(^_^2 cose)-2t;>npf (l-cos^)'. 

(20) 

The integration over momenta is governed by the positivity of the electron-positron pair 
invariant mass, P > 0. For V = 1 + 5/2 this means that we should integrate over the whole 
region —1 < cos6' < 1 for < |p'| <pc= |p|^/(4 + 5), but only over 1 > cos^ > cos6'mm — 
(-5|pp - 5|p'P + 2(1 + (5)|p||p'|)/(2|p||p'|) for pc < |p'| < |p|- The latter region gives in 
fact the major contribution to the integral. The integral for the rate of the energy loss is 
similar but with the integrand multiplied by —{E — E'). 

Performing the momentum integrals we obtain the decay rate and the rate of the energy 
loss as follows, 

Although the results are parametrically similar to the ones in Ref. [1], the numerical coef- 
ficients turn out to be model dependent. We have that for the model I, 

ai = ^(^;-l)3(^; + l)(53 + 20^;-5^;2)^^5^ (23) 

a'l = ^{v-lf{v + l)v{67 + 28v-7v^)c:^^6\ (24) 
and for the model II, 

an = ^{v' - m5v' + 19) ^6^ (25) 

a'n = ^{v'-l)M7v' + 23)^^6\ (26) 

For comparison, the results in Ref. [1], which are obtained with the standard "Lorentz 
invariant" expression for the squared matrix element, are 



{v^ - If ^ ]_ 

Uv^ ~ 14 
25 (v^ - If 25 



a' = t^^^s^ (28) 

448 t, 448 ^ ' 

We find that in all the cases, the decay rate and energy loss are proportional to S^. This 6^ 
dependence in our models reminds us of the argument in Ref. |[lj based on the kinematics 
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that a superluminal neutrino gets an "effective" mass such that the decay rate is propor- 
tional to 5^E^ . However, since the numerical coefficients are model dependent, it might be 
possible to construct models of Lorentz violation that allow for a reduction or cancellation 
of 5^ terms. 

We obtain the lifetime of a superluminal neutrino as compared to the results in [1]: 
r = is 1.76(1.25) TcG for the model 1(11). The mean fractional energy loss due to a 
single pair emission is E^^{dE / dx) /T ^ 0.81(0.78) for model I(II), which is very similar 
to 0.78, the value given in [Ij. The terminal energy of the superluminal neutrino is given 



For the OPERA baseline of 730 km, we have 13.9 GeV and 13.1 GeV for the models I and 
II, which is numerically very close to the value 12.5 GeV in |T]. 

5 Conclusions 

We calculated the decay rate and the energy loss of a superluminal neutrino in two models 
where the Lorenz violation is introduced in the kinetic term of the neutrino in the action. 
We found that due to the change in the form of the solutions of the free field (Weyl) 
equation for neutrinos, the modified spin sum rules must be used for the calculation of 
the matrix element. We also found that the final result depends explicitly on the form of 
the Lorentz violation in the action. In the analyzed models, the energy loss by electron- 
positron bremsstrahlung still makes the models incompatible with the observation of high 
energy neutrinos at OPERA (and ICARUS as well as IceCube [H [21|), but it advises us 
for a very careful calculation of the neutrino decay rate in more complicated models (for 
example, models with energy dependent velocity or modified velocities of electrons). A 
model with cancellation of 5^ contribution in the decay rate (if it exists) may evade the 
neutrino decay constraint. 
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